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Possibility of flat-band ferromagnetism in quantum dot arrays is theoretically discussed. By 
using a quantum dot as a building block, quantum dot superlattices are possible. We consider dot 
\ arrays on Lieb and kagome lattices known to exhibit flat band ferromagnetism. By performing 

£SJ ■ an exact diagonalization of the Hubbard Hamiltonian, we calculate the energy difference between 

the ferromagnetic ground state and the paramagnetic excited state, and discuss the stability of the 
ferromagnetism against the second nearest neighbor transfer. We calculate the dot-size dependence 
of the energy difference in a dot model and estimate the transition temperature of the ferromagnetic- 
paramagnetic transition which is found to be accessible within the present fabrication technology, 
^s^j , We point out advantages of semiconductor ferromagnets and suggest other interesting possibilities 

of electronic properties in quantum dot superlattices. 
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I. INTRODUCTION 



Recent progress in the fabrication technology in nano-scale ( "nano-technology" ) has enabled us to make various 
types of small devices using semiconductor quantum dots . Single electron transistor is one of the important exam- 
Ch ' pl es °f device application!!! When the charging energy of a small quantum dot (artificial atom) is larger than the 
temperature, electrons in the lead cannot transfer into the dpt due to the Coulomb blockade effect. By combining 
several single electron transistors, logic circuits are proposed.EI By coupling several quantum dots, it suggested that 
a qubit for quantum computation can be realized by controlling excess spins in coupled quantum dots.u Single electron 
transistor through coupled dots (artificial molecule) arc also jamposed, where a transition between the bonding and 
anti-binding states serves as a qubit of quantum computationD'El 

Arrays of ouantum dots have also been studied extensively. Quantum dot laser is one of the promising devices 
among them.El Recently, a self-organizing technique of quantum dotso enabled us to synthesize very small dots in 
ly-^ well-ordered lattices.u In quantum dot lasers, electrons do not transfer between dots in general and inter-dot coupling 
seems unimportant. Quantum cellular automatonEl is a very fascinating proposal of dot array device utilizing electron 
j— ( | transfer inside a cell of coupled dots and electrostatic force between the neighboring cells. Logical circuits using 
. quantum cellular automata have been also proposed£3 I— . 

Recently, the present authors have proposed a new type of device forming a superlattice of quantum dots. til If we 
consider a quantum daLas a building block and put it on a site of the lattice, we can create an artificial crystal having 
interesting properties .113 In therdpt lattice, we-ean design any type of lattice structure as we like regardless of the 
number of electrons in it. LiebEia and kagomeliil lattices are interesting examples of such-, artificial lattices, because 
they have a dispersion-less subband (flat band) in their single-particle band structures£3 It has been proven that, 
in the repulsive Hubbard model of these lattices, ferromagnetisrruappears. Interestingly, it has been shown that the 
ferromagnetic spin-wave has a finite stiffness below a Stoner gaplij and that the ferromagnetism is robuS|t-a|gainst a 
finite dispersion .113 Some other types of flat-band ferromagnetism have been proposed by several authors .E^TE] After 
their predictions, th ere h ave been several p ropo sals to realize flat-.band ferromagnetism based on real materials, such 



as carbon networks flrd a graphite ribbonjEfcl and Gac3 and Asc3 atomic wires. However, there has been no clear 
evidence of the observation of flat-band ferromagnetism, because it is difficult to form these lattices using these 



materials since a lattice distortion effect would destabilize the ferromagnetism when the flat-band is half-filled.cJ In 
real materials, the number of valence electrons is determined in such a way that the crystal structure is stable. Then, 
unrestricted material design is difficult in general. 

On the other hand, dot lattices do not have such disadvantages. One can design various types of lattice structures. 
The spatial position is fixed because an artificial atom is a rigid buried region in semiconductors. The dot lattice does 
not undergo structural deformations by electronic effects such as the Jahn- Teller distortion. Then, one can design 
lattice structures which do not exist in nature without worrying about the lattice instability. Moreover, the number 
of electrons in it can be changed in a controllable manner. By putting a gate electrode in the spatially separated 
region on top of the electron gas, it is easy to modify the electron filling by changing the gate voltage. 



1 



In this paper, we discuss the possibility of realizing flat-band ferromagnetism in quantum dot superlattices. In 
Section II, we consider two types of two dimensional (2D) dot lattices (Lieb and kagome lattices) and obtain single- 
particle band structures using a tight-binding approximation including the next nearest neighbor transfer. In Section 
III, we diagonalize a Hubbard Hamiltonian for these lattices and obtain the energy difference of the ferromagnetic 
ground state and the paramagnetic excited state. In Section IV, we describe our quantum dot model. By calculating 
the transfer and on-site Coulomb energy, we obtain the size dependence of the energy difference and discuss the stability 
of the ferromagnetic ground state. In Section V, we argue that the magnetization can be expected in 2D dot lattices in 
spite of the famous Mermin- Wagner theorem. We indicate some device applications for semiconductor ferromagnets 
and propose other possibilities of dot superlattices having interesting electrical properties like superconductivity. 
Conclusions are presented in Section VI. 



II. TIGHT-BINDING APPROXIMATION 



We consider two types of lattices. Figure [j](a) shows the Lieb lattice and Figure |l|(b) the kagome lattice. Lieb lattice 
is a bipartite lattice which consists of two sublattices. Sites belonging to one sublatticje is connected to sites belonging 
to another sublattice.E^I The kagome lattice is a line graph of the hexagonal lattice.E-3 We consider a Hubbard model 



{i,k)a i 



(1) 



where t is the transfer between a pair of the nearest neighboring sites i and j, t' is the transfer between a pair 
(i, k) of the next nearest neighboring sites i and k, U is the on-site Coulomb energy on the site i, c^.(c, CT ) is the creation 
(annihilation) operator of an electron on the site i with spin a —\ or J,, and n a — cf a Ci a . In the non- interacting case 
U — 0, we obtain the following tight-binding Hamiltonian equations of single-particle energy E for the Lieb lattice, 
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and for the kagome lattice, 
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where a is the inter-dot spacing, k = fcibi + fc2b2 is the wave vector expressed by the reciprocal lattice vectors bi 
and b2 (see insets of Fig. @), and c = (ci, 02,03) is the amplitudes of wavefunction at the site i = 1, 2, 3 in the unit 
cell shown in Fig. |l|. For t' = 0, Equations (^|) and (|3j) are simply given by 



E = 0, ±2t^Jcos 2 (k 1 a) + cos 2 (k 2 a) 

for the Lieb model, and by 

E = 2t, t(-l± ^3 + 2 cos(2fcia) + 2 cos(2fc2a) + 2 cos(2/cia — 2fc2a) 



(4) 



(5) 



for the kagome model. Figures ||(a) and (b) is the band diagrams for Lieb and kagome lattices. For t' = 0, dispersionless 
flat subbands are formed in the middle for the Lieb lattice and in the top for the kagome lattice. When the next 
nearest transfer is taken into account (t' > 0), the flat subbands are broken down, except for 2k\a = ir (from X point 
to M point) in the Lieb lattice. This is because the next nearest transfers between site 2 and 3 are cancelled between 
different unit cells. 
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III. EXACT DIAGONALIZATION 



We diagonalize the Hubbard Hamiltonian ([!]) for 2x2 unit cells of 12 dots enclosed by dashed lines in Fig. |l|. For 
the Lieb lattice, we use an anti-periodic boundary condition in the x direction and periodic boundary condition in 
the y direction in order to avoid unimportant finite-size effects of additional degeneracy at the cross point of the flat 
band and dispersive bands. For the kagome lattice, we use an ordinary periodic boundary condition. 

Figure || shows the total spin as a function of the number of electron N and U/t in the (a) Lieb and the (b) kagome 
lattices for t' — 0. For any U/t > 0, high spin states appear at the half-filling (N = 12) for—the Lieb lattice and 
19 < N < 22 for the kagome lattice. This result is consistent with the theorems proved by LiebE3 for the Lieb lattice 
and Mielkdi-3 for the kagome lattice. This result shows that the ferromagnetic state appears in a wide range of electron 
filling in the kagome lattice, whereas ferromagnetism appears only at the half-filling for the Lieb lattice. High spin 
states seen in Fig. || when the Fermi level is not at the flat band (N ^ 12 for the Lieb lattice and N < 16 for the 
kagome lattice) are due to an unimportant finite size effect due to additional degeneracy. 

Figure ^| shows the energy difference between the ground state with spin S = 2 and the lowest excited state with 
spin S = for various t's when the flat-band is half- filled, i.e., (a) N = 12 for the Lieb and (b) N — 20 for the 
kagome lattice. This energy difference can be regarded as a qualitative estimate of the transition temperature of 
ferromagnetic and paramagnetic transition in a macroscopic sample as will be discussed in Sec. IV. First we consider 
t! = 0. As Lieb has already pointed out,E3 the high spin state in the Lieb lattice is caused by the anti- ferromagnetic 
ordering at the half-filling. This ordering can be easily understood when U/t ^S> 1. In this limit at the half- filling, 
there is one localized electron per dot. Due to the second-order process of an electron with up-spin transferring to 
the neighboring dot having a down-spin electron, the effective exchange energy is given by —2t 2 /U. As the numbers 
of sites in two sublattices are different, the remaining total spin is responsible for the ferromagnetism. The high spin 
in the kagome lattice is caused by an effective exchange interaction of the thircLorder process cycling and exchanging 
two electrons with the opposite spins within the triangular lattice of three sites.E£l When U/t ^s> 1, the energy loss due 
to this ring exchange is of the order of t, whereas there is no energy loss between two electrons with the same spins. 
Actually, from Fig. |[ one can deduce AE ~ 4.3t 2 /U for the Lieb lattice and AE ~ 0.26i for the kagome lattice, 
which is consistent with the above argument. When U/t is small, the energy difference is proportional to U in both 
cases, i.e., AE ~ 0.16J7 for the Lieb lattice and AE ~ 0.05U for the kagome lattice. This is because the Coulomb 
repulsion raises the energy of the order of U when electrons on the flat band have the opposite spins in the lowest 
spin state, whereas, by the Pauli principle, it does not when electrons on the flat band have the same spins in the 
high spin state. 

As t' is increased, the energy difference becomes smaller or sometimes negative. For small U /t, the ferromagnetic 
ground state for t' = is easily broken down by very small t', because the subband is no longer flat for t' > as seen 
in Fig. [| Electron are filled in such way as to lower the total spin to gain the single-particle energy of the order of 
t' rather than as to align spins to gain the smaller exchange energy of order of U . As U/t is increased, the gain from 
the effective exchange energy AE overcomes the loss from the single-particle energy which is of order of t', and the 
ferromagnetic ground state becomes robust against t'. In the Lieb lattice, the ground state is always paramagnetic 
when t' > 0.6 in the range < U/t < 30, whereas, in the kagome lattice, the ferromagnetism disappears for t' > 0.1. 
This difference in the robustness against t' comes from the magnitude of the effective exchange energy AE as seen in 
Fig. ||. When U/t < 1, AE ~ 0.16J7 for the Lieb lattice is much larger than AE ~ 0.05J7 for the kagome lattice. For 
U/t S> 1, the anti- ferromagnetic exchange energy between the next nearest neighbor sites in the Lieb lattice slowly 
increases as (t'/t) 2 because t' 2 /U = (t' / 1) 2 (t 2 / U) , where as, in the kagome lattice, the ferromagnetic energy due to 
the ring exchange rapidly decreases as t'/t. 

IV. QUANTUM DOT MODEL 

To evaluate the transfer and on-site Coulomb energy for quantum dot arrays, we assume that electrons are confined 
in a two-dimensional confinement potential given by 

t r / \ i t-> \ i \ f — Tsm*u! 2 a 2 |cos(7ra;/a) cos(7n//a)l 2 for r < a/2 ,„s 

V(t) =^v(r-TLi), v(r) = j 16 Q for r > a/2 ' ( ) 

where m* is the effective mass of an electron, oj is the confining oscillator frequency, Hi is the position of the i-th dot, 
and a is the inter-dot spacing. This dot model is quite similar to that used for the square lattice. cB The confinement 
potential V(r) is shown in Fig. ||. Noting that [cos(irx/a) cos(7ry/a)] 2 ~ [(2r/a) 2 — l] 2 ~ 1 — 2(2r/a) 2 for r -C a/2, 
the "atomic" wavefunction localized in the potential v(r) is given in a good aproximation by 
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where d = 2yjh/m*u is the dot-diameter. The transfer and on-site Coulomb energies are given by 

t(R h Rj)= - y dr ( f)(R l ){-h 2 X7 2 /2m* + y(r)}0(r - Rj), (8) 



where e is the dielectric constant. Note that the minus sign of the integral in Eq. (|J) comes from the definition of 
the transfer energy in Eq. (|l|). Here, the nearest neighbor transfer t is calculated for Rj = (0,0) and Rj = (a, 0) 
and the next nearest neighbor transfer if for Rj = (a, 0) and Rj = (0, a) in the Lieb lattice and by Rj = (a, 0) 
and R 3 = (a/2, >/3a/2) in the kagome lattice. The Coulomb energy can be analytically integrated, but the transfer 
integral is evaluated numerically from Eqs. (0) and (0). Figure ^ shows the calculated transfer and the on-site Coulomb 
energy as a function of the dot-diameter. Here, we adapt the effective atomic units h 2 /2m* — e 2 /Aire = 1. In these 
units, energy and length are scaled in units of the effective Rydberg constant Ry* = 13.6 cV x (m* /m n )/(e/eo) 2 and 
the effective Bohr radius a* B — 0.53 A x (e / eq) / (m* / mo) , where mo and £o is the mass of a bare electron and the 
dielectric constant in vacuum. In our dot model, the contribution from terms like (0(O)|u(Rj)|</>(Rj)) (Rj ^ 0, Rj) 
is negligible in Eq. (|^), and the transfer energy is determined by the distance between the nearest neighboring dots, 
i.e., t ~ (</>(0)|{p 2 /2m* + v(0) + «(Ri)}|</>(Rj)). Therefore, the transfer energies for the Lieb and kagome lattices 
are almost identical within the width of drawn lines in Fig. ^. On the other hand, t' for the Lieb lattice is larger 
than that for the kagome lattice, because of the difference in the distance of the next nearest neighbor dots, i.e., 
|Rj — Rj| = v2a for the Lieb lattice and |R — Rj| = V3a for the kagome lattice. For d/a < 0.5, t' is much smaller 
than t and negligible. It is noted that t is always smaller than U when the inter-dot spacing a > 0.1, which is usually 
realized in the present fabrication technology. 

In realistic dot arrays formed by a negatively biased gate electrode depleting the underneath two-dimensional 
electron gas, the inter-dot spacing is usually fixed and cannot be changed. By modifying the gate voltage, the dot- 
diameter can be changed. To simulate this, we evaluate t, t' , and U from Eqs. (g) and (^|) for a fixed inter-dot spacing 
(a=0.5, 1, 5, 10) and calculate the energy difference as a function of dot diameter as shown by the closed symbols in 
Figure ^. The energy difference is appreciable for 0.3 < d/a < 0.7. For a > 0.5, t is always smaller than U. In this 
strongly correlated region, the energy difference AE has the dependence of AE ~ t 2 /U ~ 1/a 3 for the Lieb lattice, 
and AE ~ t ~ 1/a 2 for the kagome lattice. For d/a > 0.6, AE rapidly decreases, since if becomes significantly large 
and the ground state becomes paramagnetic. In Fig. fj], AE for if — is also plotted. When d/a > 0.6, AE for if = 
monotonically decreases but does not become negative, since the ground state is always ferromagnetic for t' = 0. It is 
noted that, when AE takes a peak around d/a = 0.5 — 0.6, the effect of t' is negligible and AE is not affected by t 1 . 

In Table §, we estimate AE for dot arrays of various sizes. We consider GaAs, InAs, and Si dots and assume a — 2d. 
For dots with a spacing of lOOnm which is available within the present fabrication technology, AE is of the order of 
several hundreds mK and we can expect that ferromagnetism can be observable in the dilution temperature region. 
For dots of spacing of 5nm, AE is as high as a few tens of K. 



V. DISCUSSION 



There are several advantages of using semiconductors in making artificial crystals. First, the lattice structure can 
be widely chosen. One can fabricate a lattice structure which does not exist in nature. Second, inter-dot coupling 
and the electron filling can be separately modified. This is possible if the inter-dot coupling is modified mainly by 
the front-gate electrode on top of the two-dimensional electron gas and the electron filling is modified mainly by the 
back-gate electrode. The controllability of electron filling enables us to switch ferromagnetism on and off. It would 
be better to use the kagome lattice in order to switch the ferromagnetism because the ferromagnetic ground state 
appears in a wide range of electron filling as shown in Fig. ||(b). The effect of a magnetic field in the kagome lattice 
is also very interesting because the flat band is destroyed by the threaded magnetic flux. It has been found that a 
giant negative magpetoresistance and ferromagnetic-paramagnetic transition induced by a magnetic field occur in the 
kagome dot latticeO Recently, the present, authors have proposed an simple way to realize the kagome dot lattice 
within the present fabrication technology.t3 It has been shown that a network of quantum wires effectively acts as 
a kagome lattice where electrons are well localized at the cross points of two wires. This kind of kagome network 
has been already available in quatum wires formed by a selective area growth techniques This method of making a 
kagome dot lattice provides us a chance to observe the flat-band ferromagnetism in experiment. 
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One may think that there should be no long-range -order such as ferromagnetism in two dimensions at finite 
temperatures according to the Mermin- Wagner theorem.E3 However, it has been shown that the spin-spin correlation 
length £ in the spin 1/2 Heisenberg model on a 2D square lattice exponentially grows like £ ~ exp(J/T) as decreasing 
temperature although the (anti- ferromagnetic) spin-spin correlation decays like (So ■ S r ) oc exp(— r/£)H3 Then, the 
spin-spin correlation length £ can be macroscopically large at low temperatures. For example, at T — 0.1J, £ is 
10 4 times larger than the lattice spacing. Moreover, the the spin-spin correlations in the 2D XY model d$£3y very 
slowly with distance r like r -1 / 4 below the Kosterlitz-Thouless-Berezinskii (KTB) transition temperature,EjEa and a 
measurable magnetization is always present in any realizable system having a finite size. Actually, below the KTB 
transition temperature Tkt, the thermal-averaged magnetization per site is shown to be M(N,T) — (2N)~ T / S7rJ , 
where J is the coupling constant and N is the number of sites.B At T = T KT = ttJ/2, M (N, T K t) = (27V)- 1 / 16 - 0.54 
for N — 10 4 |-^rhe effective Curie temperature for a finite-size system where the finite magnetization appears is similar 
to Tkt ~ </■" These two results suggest that a finite magnetization can be expected even in two-dimensions as long 
as the sample size is smaller than the spin-spin correlation length. Moreover, considering J ~ AE, our estimation of 
the transition temperature in Fig. [?] is qualitatively justified. 

Although a finite magnetization is expected in two dimensions, its value will be quite small in 2D dot arrays, 
since the lattice constant (> 10 nm) of dot arrays is about more than a hundred times larger than that of the 
conventional ferromagnetic materials and the expected magnetization per area will be significantly reduced. To 
measure the magnetic moment directly, a very sensitive detector such as a SQUID magnetometer would be required. 
A more convenient way to detect the magnetization directly will be to measure the anomalous Hall resistivity which 
is proportional to the magnetization of the sample and is added to the normal Hall resistivity proportional to the 
external magnetic field. On the other hand, magnetoresistance measurement is rather an indirect way to detect 
the magnetization. When the magnetic field is increasedj the insulating ferromagnetic state turns into the metallic 
paramagnetic one in the kagome lattice at the half-filling.E3 The magnitude of the magnetization could not be estimated 
only from the magnetoresistance measurement, although the ferromagnetic-paramagnetic phase transition could be 
detected. 

The advantage of using semiconductors in realizing ferromagnetism exists not only in making semiconducting 
"permanent magnets" . As we mentioned, the controllability of the magnetic property by changing the electron filling 
and the magnetic field will make dot lattices useful in electronic devices such as memories, sensors, and magnetic 
heads since only semiconductor materials such as Si and GaAs are contained. Magnetic devices can be fabricated 
without using any magnetic elements such as iron and manganese which are incompatible with the conventional LSI 
fabrication process. 

By extending the idea of dot superlattices, one can think of other interesting possibilities of artificial materials. 
Thanks to the rapid progress in the semiconductor nano-technology, we can expect that various interesting electric 
properties which have been observed in conventional materials may be realized also in dot superlattices. One of the 
most fascinating examples is high-temperature (high-T c )-superconductivity. High-temperature <i-wave superconduc- 
tivity in a repulsive Hubbard model has been predicted.EZI Although the energy scales or the transition temperature 
is a hundred times smaller than the conventional Cu02 high-T c superconductors (T c ~ 100 K), superconductivity 
in semiconductor dot arrays might be possible, because the estimated transition temperature using the predicted 
expression is T c ~ O.Olt ~ 1 K for a = 2d = 10 wa-GaAs dots. It has been shown that, in other types of lattices, 
the transition temperature becomes much higher Ptt 1 It would be very interesting if superconductivity (or at least 
the KBT transition) could be realized in semiconductors. Other types of lattice structures would also be fascinating, 
such as a ladder structure realized in coppeiijaxide materials, where various interesting properties have been observed 
such as the spin gap and superconductivityOEj Optical properties of the dot arrays would also be interesting, since 
the large density of states in the flat band will significantly affect photo-luminescence or laser characteristics. 

VI. CONCLUSIONS 

Flat-band ferromagnetism in quantum dot arrays has been theoretically discussed. We considered dot arrays on 
the Lieb and kagome lattices which are known to exhibit flat-band ferromagnetism. The tight-binding calculation 
showed that the next nearest neighbor transfer t' destroys the flat subband. We performed the exact diagonalization 
of the Hubbard Hamiltonian and calculated the energy difference between high-spin and low-spin states. This energy 
difference represents a qualitative estimate of the transition temperature of ferromagnetic and paramagnetic states in 
a macroscopic sample. The energy difference becomes smaller as the next nearest neighbor transfer t' increases. It 
was shown that, although the ferromagnetic ground state is easily broken down by t' in the weak correlation region 
(t ^> U), it is robust against t' in the strong correlation region (U ^> t). 

We calculated the size dependence of the energy difference in a realistic dot model. We found that, although the 
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next nearest neighbor transfer destroys the ferromagnetism when the dot diameter approaches the inter-dot spacing, it 
does not affect the peak value of the energy difference or the transition temperature when the dot diameter decreases. 
We argued that the flat-band ferromagnetism can be observable in dot arrays fabricated using the present technology. 
We suggested other interesting possibilities of artificial material design using quantum dot superlattices. 
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FIG. 1. Lieb lattice (a) and kagome lattice (b). Solid lines represent the nearest neighbor transfer t. Dotted and dashed 
lines indicate one and 2x2 unit cell(s) respectively. A unit cell contains three sites i = 1, 2, and 3. 
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FIG. 2. Single-particle band energy E for (a) Lieb lattice and (b) kagome lattice for the next nearest neighbor transfer t' = 
(solid line) and t' > (dashed line) calculated using the tight-binding approximation, where E and t' are normalized in units 
of the nearest neighbor transfer t. Insets: Brilouin zone, bi and i>2 are the reciprocal lattice vectors. 
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FIG. 3. Total spin as a function of the number of electrons < N < 24 for U/t = 10 3 , 5, 10, 15, and 20 for (a) the Lieb 
lattice and for (b) the kagome lattices when t' = 0. 
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FIG. 4. Energy difference E(S = 0) and E(S — 2) as a function of U/t in (a) the Lieb (N = 12) and (b) the kagome lattices 
(N — 20) for various values of t'. 
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FIG. 5. Potential profile of the quantum dot model in Eq. (j6|) for (a) the Lieb and (b) the kagome lattices. Length is 
normalized in units of the inter-dot spacing a. 
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FIG. 6. The nearest and the next nearest neighbor transfer energies t, t' and on-site Coulomb energy U in the dot model as 
a function of d/a. Energy and length are scaled in units of the effective Rydberg constant Ry* and the effective Bohr radius 
a* B (see text). 
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FIG. 7. Energy difference AE = E(S = 0) — = 2) as a function of dot-diameter d for (a) the Lieb and (b) the kagome 
lattices for various inter-dot spacing a in the dot model represented by solid symbols. Open symbols represent the energy 
difference when we neglect the next nearest neighbor transfer t' . 
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TABLE I. Estimated energy difference AE — E(S = 0) — E(S = 2) (kelvin) between the ground state with spin 2 and the 
lowest excited state with spin for dot arrays of the kagome lattice of typical inter-dot spacing a =5, 10, 50, lOOnm and dot 
diameter a/2. We consider GaAs dots (m* /mo = 0.067, e/eo = 12.4, Ry* = 6meV, a* B = lOnm), InAs dots (m* /mo = 0.02, 
e/e = 12.4, Ry* = 1.8meV, a* B = 34nm, and Si dots (m* /mo = 0.2, e/eo = 12, Ry* = 19meV, a* B = 3nm). 
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AE (kelvin) for Lieb lattice 
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